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ABSTRACT 

The dwarf spheroidal (dSph) satellites of the Local Group have long been thought to 
be simple spheroids of stars embedded within extended dark matter halos. Recently, 
however, evidence for the presence of spatially and kinematically distinct stellar pop- 
ulations has been accumulating. Here, we examine the influence of such components 
on dynamical models of dwarf galaxies embedded in cold dark matter halos. We begin 
by constructing a model of Andromeda II, a dSph satellite of M31 which shows evi- 
' dence for spatially distinct stellar components. We find that the two-component model 

\Q , predicts an overall velocity dispersion profile that remains approximately constant at 

OO ' ^ 10 — llkms -1 out to ~ 1 kpc from the center; this is despite wide kinematic and 

C^> spatial differences between the two individual components. This prediction can be 

^sO . validated by detailed spectroscopic analysis of this galaxy. The presence of two com- 

ponents may also help to explain oddities in the velocity dispersion profiles of other 
^ (— i dSphs; we show that velocity dispersion profiles which appear to rise from the center 

O outwards before leveling off — such as those of Leo I, Draco, and Fornax — can result 

I ' from the gradual transition from a dynamically cold, concentrated component to a 

second, hotter, and more spatially extended one, both in equilibrium within the same 
dark halo. Dwarf galaxies with two stellar components generally have a leptokurtic 
line-of-sight velocity distribution which is well described by a double Maxwellian. This 
may be contrasted with other dynamical explanations such as a radially-dependent 
anisotropy in the stars' orbits. Interestingly, we find that multiple equilibrium compo- 
nents could also provide a potential alternative origin for "extra-tidal" stars (normally 
ascribed to tidal effects) in situations where corroborating evidence for tides — such as 
elongation of the main body of the dwarf in the orbital direction or velocity gradients 
across its face driven by protruding tidal tails — may be lacking. 
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1 INTRODUCTION 

Dwarf spheroidal (dSph) galaxies are extremely faint sys- 
tems whose relatively large size and velocity dispersion com- 
bine to make them premier examples of systems likely dom- 
inated by dark matter. They contain little or no gas, and 
individual stars are the only available tracer of their struc- 
ture and kinematics. Although long regarded as relatively 
simple, single-component spheroids of stars, this view has 
been challenged by the advent of panoramic imaging tech- 
niques and multi-object spectrographs in large telescopes. 

The analysis of the color-magnitude diagrams (CMDs) 
of dSphs shows, in many systems, compelling evidence for 
multiple and protracted episodes of star formation (eg. 
Grcbcl 1997; Mateo 1998, and references therein), as well 
as for the presence of spatially distinct stellar populations 



(|Da Costa et al.lll996l : lHarbeck et alj|200ll ). These popula- 
tions often spill beyond the nominal "tidal radius" estimated 
from K ing-model fits to t he inner surface brightness pro- 
file (eg. iMufioz et al.ll2006l and references therein), a result 
which has brought into focus questions regarding the role 
of Galactic tides in limiting the spatial extent of dSphs, as 
well as the m ass and spatial extent of their surrounding dar k 
matter halos l| Johnston et al.|[20o3 : iPenarrubia et al.ll2007l '). 

Spectroscopic studies have also uncovered remark- 
able complexity in the dynamics of dSphs. As suf- 
ficient data become available, a complex picture has 
arisen where kinem atic oddities such as cold clumps 
l|Klevna et alj 20031). chemodynamically distinct stellar 
populations (iTolstov et al.l 12004 iBattaglia et~all 120061 : 



llbata et al l l2006h, and hints of ri sing dispersion profiles 
( Wilkinson et al.ll2004l ; IMufioz et al]|2005l : iKoch et al.ll2007l : 
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I Walker et all 120061 ), are common. This is at odds with the 
"natural" expectation for well-mixed stellar systems, where 
the velocity distribution is expected to be a smooth func- 
tion of radius. Observations thus suggest that it may be 
time to relax the standard assumption that dSphs are well 
described by a single well-mixed stellar system embedded 
within a dominant dark matter halo. 

In IPenarrubia et al.l (|2007t ) . we examined the observa- 
tional properties of single component stellar systems (de- 
scribed by King models) embedded in massive cold dark 
matter CDM) haloes (described by Navarro, Frenk & White 
(1996, 1997, hereafter NFW) profiles), and applied these 
models to the known population of Local Group dSphs. 
We examine here the general effect of allowing for multi- 
ple stellar components in these dynamical models. We begin 
in §2 by motivating the modeling technique using the case 
of Andromeda II (And II; a dSph companion of M31), and 
then examine in general the velocity dispersion and surface 
density profiles of two-component stellar systems embedded 
within a CDM halo. We conclude in §3 with a discussion of 
our results and their wider implications. 



2 DYNAMICS OF MULTIPLE STELLAR 
COMPONENTS 

2.1 The case of Andromeda II 

And II is a dSph satellite of M31 lo cated at a distance 
of ~ 650 k p c from the Milky W a y |McConnachie et al.l 
12004 120051 '). iMcConnachie fc Irwinl $2006) present its sur- 
face brightness profile based on Isaac Newton Telescope wide 
field imaging, and suggest that a discontinuity in the pro- 
file at r ~ 2' may be evidence that it possesses a secondary 
"core" component. 

Deeper, follow-up im aging of And II was obtained by 
IMcConnachie et al.l l|2007h using the Subaru Suprime-Cam 
wide field camera. The black squares in Figure 1 reproduce 
the overall radial profile of And II derived using these data. 
The deeper data allow the derivation of the radial profiles of 
distinct stellar populations; the bottom two profiles shown 
in the same panel trace the spatial distribution of horizontal 
branch stars (HB) and the reddest red-giant branch (RRGB) 
stars. Clearly, the HB stars are more spatially extended than 
the RRGB stars, which trace an exponentia l profile. As dis- 
cussed in detail in lMcConnachie et al.l ([2007), it is very likely 
that these two profiles represent the density distribution of 
two distinct stellar components in And II. Not only does 
an appropriately weighted sum of these two profiles fits the 
overall number density profile remarkably well, as indicated 
in Figure 1, but it also explains significant differences seen in 
the stellar po pulations at small and larg e radii. The reader 
is referred to IMcConnachie et al.l (I2007T ) for more detail on 
these data and on the construction of these profiles. 

Large differences in spatial distributions of stars should 
manifest themselves kinematically, and this can in principle 
be used to constrain the properties o f the overall poten- 
tial. As described by IPenarrubia et alj l|2007h . one expects 
a direct link between the velocity dispersion and the spatial 
distribution of stars orbiting within a CDM halo. A stellar 
component deeply embedded within the central regions of a 
CDM halo is expected to have a low central velocity disper- 
sion and a dispersion profile that remains flat well outside 



its characteristic core radius. More extended components, 
on the other hand, should have higher central velocity dis- 
persions and more steeply declining dispersion profiles. 

Detailed k inematic st u dies o f And II are currently lack- 
ing. However, ICote et"al1 (|l999h measure a central veloc- 
ity dispersion of 9.3 ± 3 kms - based on 7 stars. Given the 
expected dominance of the central component in the inner 
few hundred parsecs of And II, it is likely that this mea- 
surement is representative of the centra l velocity dispersion 
of the exponential component fou nd by McConnachie et al.l 
i|2007l ). Following the technique of IPenarrubia et all l|2007h . 
we adopt this value and the appropriate value for the 'core' 
radius of the central component (r c ~ 370 pc) to deter- 
mine the properties of a ACDM halo at z — compatible 
with these constraints. The resulting NFW profile may be 
fully characterized by the peak of the circular velocity curve, 
(V m ax, J"max): we find (Vmax = 17.4 km s" 1 , r max = 3.06 kpc). 

The velocity dispersion profiles of both components can 
then be computed, assuming isotropy, by solving Jeans' 
equations separately for each component within the NFW 
potential. The predicted velocity dispersion profiles for 
And II are shown in the top right panel of Figure 1. As ex- 
pected, the exponential component has a lower central veloc- 
ity dispersion than the extended component and its profile 
declines by only ~ 1 kms" 1 out to R ~ 1 kpc. The more ex- 
tended component, on the other hand, has a steeply declin- 
ing profile, dropping from a central value of ~ 14.3 km s -1 
to ~ 7 kms -1 at R ~ 2 kpc. For comparison, the circular 
velocity profile of the NFW halo is plotted as a gray dashed 
line as a function of actual, not projected, radius. 

The overall velocity dispersion profile of And II is given 
by the weighted sum of the profiles of the components. It 
has a central velocity dispersion of ~ 11 kms" 1 which rises 
slightly and then declines slowly; at R ~ 1 kpc, the velocity 
dispersion is ~ 10 kms -1 , and at R ~ 2 kpc the dispersion 
has dropped to ~ 7 kms -1 . Varying the adopted values of 
(Vmax,r ma x) by using t he different constraints discussed in 
IPenarrubia et ail l|2007h does not alter the central velocity 
dispersions significantly. The general shape of the overall 
profile of And II is similar to that expected for a one com- 
ponent galaxy, but here it is a result of the varying radial 
contribution of two dynamical components; this is very dis- 
tinct from the standard interpretation of a flat dispersion 
profile arising from only one stellar component 

Our predictions for And II are verifiable using current 
generations of multi-object spectrographs by obtaining spec- 
troscopic metalliciti e s of r ed giant branch stars in And II. 
IMcConnachie et al.l l|2007l ) conclude that the two compo- 
nents of And II have distinct metallicities, and so it should 
be possible to divide the red giant branch stars based on 
their met allicity. This techniq ue has been successfully e m- 
ployed bv lTolstov et al.l (|2004T ) and lBattaglia et all \200dh to 
uncover kinematically distinct components in Sculptor and 
Fornax, respectively. Therefore, the model which we propose 
for And II is eminently falsifiable. 

One may also resort to the distribution of line-of-sight 
velocities. If And II is largely made up of two independent 
equilibrium components of widely different dispersion, the 
velocity distribution will approach a double Maxwellian dis- 
tribution whose kurtosis (/i4//i| — 3) will be positive (lep- 
tokurtic) everywhere where both components contribute in 
comparable amounts to the density profile. This is shown in 
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Figure 1. Left panel: the density profile of Andromeda II derived by McConnachie, Arimoto & Irwin (2007) (black squares). This is 
consistent with the sum of two spatially distinct stellar components; an extended component traced by horizontal branch stars (HB, blue 
diamonds/ dotted line) and an exponential component traced by the reddest red-giant branch stars (RRGB, red triangles/dot-dashed 
line). Upper right panel: The predicted projected velocity dispersion profile for Andromeda II (solid curve), together with the velocity 
dispersion profiles of the individual components, shown by dot-dashed and dotted curves. Lower right panel: the kurtosis of the line-of- 
sight velocity distribution as a function of radius for this model. Dwarfs with multiple equilibrium stellar components are expected to 
have a leptokurtic velocity distribution at all radii where both components contribute significantly. 
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Figure 2. Left panel: Projected density profile of a system consisting of the superposition of two King models with mass ratios as 
labeled. Top right panel: Projected velocity dispersion profiles of the two stellar components and the overall projected velocity dispersion 
profile (line styles are as in the bottom right panel). Bottom right panel: The relative contribution of the more extended component to 
the total projected stellar density, shown as a function of radius for various choices of the mass ratio M2/M1. 



the bottom right panel of Figure [TJ where the kurtosis peaks 
at the center and declines gradually in the outer regions as 
the contribution from the exponential component drops. The 
uncertainty in the kurtosis measurement is where 
TV is the number of stars in the sample. Thus to measure 
the predicted value of the kurtosis in the central regions of 
And II to ~ 3 — a will require of order 200 stars; large kine- 
matic sample are crucial for accurate kurtosis measurements 
in dwarf galaxies. 



2.2 Two-component systems 

As a more general illustration of the interplay between cr p (R) 
and multiple components, we show in the left panels of Fig- 
ure 2 the density profile corresponding to two King models 
embedded within an NFW halo. We use King models to de- 
scribe the stellar distributions here because they are a con- 
venient means of parameterising the density profile and have 
been used extensively is the literature; no additional physics 
is implied by their usage. We set the core and tidal radii of 
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the more extended component to r ci = 400 pc, r tl = 2000 pc 
(dotted lines), and r C2 = 100 pc, r t2 = 400 pc (dot-dashed 
lines) for the more compact one. The mass ratio of the two 
components varies from M2/M1 = 0.1 in the top panel to 
M2/M1 = 10 in the bottom panel. The overall stellar density 
distribution is shown as a solid line in each panel; that of 
the top panel is unlikely to be observationally distinguished 
from a single component system in the absence of other data; 
the one in the bottom panel shows a "bump" in the outer 
density profile, not unlike those reported in som e dSphs (eg. 
Ilrwin fc Hatzidimitriou1ll995l ; ISohn et alj|2006h . 

The top right panel of Figure 2 shows the line-of-sight 
velocity dispersion profile, a p (R) , corresponding to these 
two components (embedded in an NFW halo with r ma x = 2 
kpc and V ma x = 30 km s _1 ). The bottom right panel shows 
how the importance of the hotter, more extended component 
(#1) varies with R for various choices of M2/M1. Note that, 
as in the case of And II, the more extended component (#1) 
has a declining velocity dispersion profile, in contrast with 
that of the more concentrated one (#2), which remains more 
or less flat all the way to its tidal radius, r t2 . 

The detailed shape of the overall velocity dispersion pro- 
file depends on the relative contribution of the two compo- 
nents and on their velocity difference. The various lines in 
the top panel of Figure [2] illustrate this for various values 
of the ratio M2/M1 (see key in the bottom right panel). 
When one component dominates everywhere the velocity 
dispersion profile is indistinguishable from a single compo- 
nent model. As variations in the relative contribution of the 
two components become more significant, however, the ve- 
locity dispersion profile starts to deviate from that for a 
single component system. The exact shape of the resulting 
<j p (R) profile depends closely on Ei/ (Ei + E2); o p (R) may 
therefore set additional constraints on the spatial properties 
in two component models. 



3 DISCUSSION 

In two component models, the centrally concentrated com- 
ponent is dynamically colder and has a decreasing influ- 
ence on the overall dispersion profile as a function of ra- 
dius. In the extreme case where each component dominates 
at different radii, we transition from a cold to a hot ve- 
locity dispersion profile. It is intriguing to note that such 
features in the velocity dispersion profile may not be un- 
usual in Local Group dSphs. This is illustrated in Figure [3j 
where we show data for Leo I, Fornax, and Draco, com- 
piled from the literature. In each of these three cases, a p 
appears to rises by about a factor of ~ 2 from its central 
value before starting to de cline again. In t he ca se of For- 
nax, we refer the reader to iBattaglia et al.l (|2006h . The ve- 
loc ity dispersion profi l e of D raco has previously been fitted 
by iMashchenko et al.1 (2006) using a single stellar compo- 
nent with a radially v arying orbita l aniso tropy embedded 
m a massive halo, and iKoch et al.l (|2007h imply a similar 
anisotropy profile is required for Leo I if it resides in an 

NF W halo. 

iKlimentowski et al.l (|2006T ) suggest that that these 
bumps in the velocity dispersion profiles are due to tidal ef- 
fects. However, we consider this un likel y for two reasons (i) 
as concluded bv lRead et al.l (|2006al ) and IKlimentowski et al.l 
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Figure 3. Velocity dispersion profiles of three Local Group dSphs 
which show a rise in their velocity dispersions from the center 
outwards before leveling off, taken from the literature. Possible, 
but degenerate, decompositions of these profiles into a cold and 
hot component are shown to illustrate the general plausibility of 
two dynamical components instead of one. 



(2006), tides do not appear to affect dSph satellites inside 
of 1 kpc, and so are unlikely to produce observable fea- 
tures at radii of typically 100 — 400 pc (ii) the presence 
of tidally stripped stars manifest themselves as monotoni- 
cally incr easing features in the velocity dispersion p rofiles 
of dSphs l|Read et al.ll2006al ; IKlimentowski et al.ll2006l '). The 
features highlighted in Figure 3, however, rise and then level 
off or decline. Thus we conclude their origin is unlikely tidal. 

For each galaxy in Figure 3 we have shown a possible 
decomposition of the velocity dispersion profile into a cold 
and hot component. These curves are inevitably degenerate, 
given that many parameters need to be fitted to describe 
two stellar components and a dark matter halo. In addition, 
we have not attempted to determine the uncertainties due 
to observation errors or bin sizes since these profiles have 
been taken direct from the literature. Some caution is clearly 
required; for example, recent data for Leo I do not appear 
to show such as pronounced rise in the dispersion profile as 
found by Koch et al. (M. Mateo, private communication). 
As such, these curves merely demonstrate the plausibility of 
multiple components in these systems and illustrate a few 
fairly general conclusions that can be drawn. 

The first is that the rise in a p from the center out- 
wards reflects the rise in the relative contribution of the 
more extended component, a fact that may be used to con- 
strain the relative size of the two components and to pin- 
point the radius where both components contribute more 
or less equally to a p (R). This is where deviations from a 
simple Maxwellian-like velocity distribution are likely to be 
maximized. At this position, the line-of sight velocity distri- 
bution would be best characterized by a double Maxwellian 
distribution, distinguishing it from other proposals such as 
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a radi ally- varying orbital anisotropy fsee lMashchenko et alj 
(2006) for application of this idea to Draco). 

The ratio of central-to-maximum velocity dispersion 
provides a lower bound to the ratio of the central veloc- 
ity dispersions of the two components. Interpreting Leo I as 
a two-component King model, the data shown in Figure 
suggest that the central component has a central (projected) 
velocity dispersion of ~ 6-7 km s -1 , and that the outer com- 
ponent has cr P2 (0) ~ 13-14 km s _1 . If Leo I is a multiple 
component system, we expect that the line-of-sight veloc- 
ity distribution of stars at about R ~ 150 — 200 pc should 
reveal evidence for the double Maxwellian distribution dis- 
cussed previously. Similar testable inferences can be made 
for Fornax and Draco. 

We note that a multiple component system is not re- 
quired to exhibit differences in the stellar populations of 
its components. For example, one could build a dwarf by 
merging two progenitors of very different spatial extent and 
kinematics but similar age and metallicity. The remnant of 
such a merger would appear homogeneous to a CMD analy- 
sis but would retain the signature of the two progenitors in 
phase space. Within this context, evidence for a hierarchical 
origin of a dSph might be best appreciated in the dynamical 
evidence rather than in the CMDs. 

Perhaps the most common deviation from single com- 
ponent models is the presence of stars that are clearly associ- 
ated with the dSph, but which lie beyond the putative limit- 
ing radius of a King model fit. The presence of these "extra- 
tidal" stars can be caused by dynamical heating in the pres- 
ence of the tidal field of the Galaxy and has attracted con- 
siderable attention because of its potential to constrain the 



mass and extent of a dS ph 's dark halo (e.g., Johnston et alj 
120021 ; iRead et aij|2006al ltt iKlimentowski et al.ll2006h . 

As illustrated in Figure 2, the presence of a low density, 
extended second component, could in some instances be mis- 
taken for "extra-tidal" stars, although in this scenario the 
stars are actually bound and in equilibrium with the dark 
matter halo. Note also that, in this interpretation, we would 
expect a v to decline in this region, since these stars belong to 
the extended component. Thus, a multiple component sce- 
nario offers an alternative explanation for "extra-tidal" stars 
in situations where corroborating evidence for tides — such 
as elongation of the main body of the dwarf in the orbital 
direction or velocity gradients across the face of the dwarf 
driven by protruding tidal tails — may be lacking. 

If the multiple component scenario is correct, a natural 
question which arises from this discussion is: how did these 
systems develop and preserve such complex structures! A 
spatially varying star formation history, where subsequent 
epochs of star formation occurred in different volumes tha n 
previous ones, has been sugge s ted bylKawata et al.l (2006). 
Alternatively, iBattaglia et alj (|2006l ) suggest that there is 
tentative evidence for non-equilibrium kinematics in their 
data for Fornax, which might imply a merger-driven sce- 
nario. However, it is currently too soon to distinguish be- 
tween these and other scenarios. Given the observational dis- 
covery of multiple structural components in dSph galaxies, 
it is important that the consequences for dynamical models 
of these systems are fully explored. 
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